Strong-field approximation for intense-laser atom processes: the choice of gauge 
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The strong-field approximation can be and has been applied in both length gauge and velocity 
gauge with quantitatively conflicting answers. For ionization of negative ions with a ground state 
of odd parity, the predictions of the two gauges differ qualitatively: in the envelope of the angular- 
resolved energy spectrum, dips in one gauge correspond to humps in the other. We show that the 
length-gauge SFA matches the exact numerical solution of the time-dependent Schrodinger equation. 
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Quantum mechanics is gauge invariant: it is easily 
proven that a given physical quantity can be evaluated 
in any gauge with the same result |]J. In nonrelativis- 
tic quantum mechanics, when the dipole approximation 
is adopted, the interaction of an atom with a time- 
dependent field such as a laser field is usually described 
in either one of two gauges: the length gauge (L gauge) or 
the velocity gauge (V gauge). In numerical solutions of 
the time-dependent Schrodinger equation (TDSE), gauge 
invariance has been confirmed many times. In analyti- 
cal work, however, some approximations almost always 
have to be adopted. There is no formal reason of why af- 
ter such approximations the resulting theory should still 
be gauge invariant. Indeed, the lack of gauge invariance 
after what seems like very well justified approximations 
has driven many a researcher to despair [2| ■ 

In this paper, we will address one of the most glar- 
ing manifestations of this "gauge problem" : the lack of 
gauge invariance of the strong-field approximation (SFA) 
in intense-laser-atom physics 0] ■ The SFA underlies al- 
most any analytical approach to total ionization rates, 
above-threshold ionization, high-order harmonic genera- 
tion, and nonsequential double ionization, both of atoms 
and of molecules. Briefly, it assumes that the initial 
bound state of the atom or molecule is unaffected by the 
laser field while the final state, which is in the contin- 
uum, does not feel the presence of the binding potential. 
The lack of gauge invariance of the SFA has been noted 
many times; see, e.g., Ref. |4|. Comparisons that have 
been carried out indeed have exhibited significant dis- 
agreements between the results obtained from L gauge 
and V gauge Different authors have preferred dif- 
ferent gauges. The question of which gauge is superior 
for which problem has often been raised, but never led 
to any consensus about its answer. Below, we will give 
an answer for the case of a short-range binding poten- 
tial, where the SFA is expected to be most accurate [6|, 
by comparing the SFA in L gauge and V gauge with the 
numerical solution of the TDSE. 

For a fixed nucleus and in the single-active-electron ap- 
proximation, where the effects of all electrons but one are 
absorbed into an effective binding potential, the complete 
Hamiltonian in the presence of an external electromag- 



netic field can be decomposed as 

H x (t) = H Q + H Ix {t), (1) 
where the subscript x specifies the gauge (x — L, V) and 

P 2 

H = — + V(r) with p = -iV. (2) 
2m 

This operator contains the binding potential V(r) and is 
independent of the gauge. With the dipole approxima- 
tion, which neglects the space dependence of the electric 
field and the vector potential, so that E(r, t) — ► E(t) 
and A(r, t) — > A(i), respectively, the electron-field inter- 
action operator has the following forms in length gauge 
and velocity gauge: 

A free electron (no binding potential) in the presence of 
the laser field is governed by the Hamiltonian 

P 2 

H Fx (t) = ^-+Hi x (t). (4) 
Zm 

The time-evolution operator of the total Hamiltonian JJJ 
satisfies the Dyson equation 

t 

U x (t,t') = U {t,t')-i J dTU x (t,T)H Ix ( T )U (T,t'), (5) 
t' 

where Uo(t, t') denotes the time-evolution operator of the 
Hamiltonian J3J). 

The exact ionization amplitude from an initial bound 
state \ipo(t)) = 0)exp(«/p£) with ionization potential I p 
to a final continuum state \ip P {t)), both defined by the 
Hamiltonian Hq, is 

M p = lim {i; p (t)\U x (t,t')\Mt'))- (6) 

t — >OQ,t' — > — OO 

We assume that the laser field be turned off in the limits 
of t — ► oo and t' — > — oo and that A(oo) = A(— oo) = 0. 
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Gauge invariance then implies that M p be gauge invari- 
ant, and indeed this can easily be verified explicitly. The 
SFA is obtained if we insert the Dyson equation © into 
the ionization amplitude ©. The first term, which comes 
from Uo(t, t'), cancels since the initial and the final state 
are orthogonal, and we are left with Q 



with the action 



M D = — i lim 

v t—>oo 



dT(iP p (t)\U x (t,T)H lx (T)\ip (r)), (7) 



which is still exact. 

In the argument that follows we restrict ourselves for 
the sake of transparence and simplicity to "direct" elec- 
trons, i.e., those that after the initial ionization never 
again feel the binding potential. In order to obtain 
the transition amplitude for the direct electrons, we re- 
place in Eq. Q the exact state at time r, which is 

(ip p (t)\U x {t, r), by the Volkov state (^> P x (t)| [given be- 
low in Eq. ((5J)] where the interaction with the binding 
potential is neglected. This yields the well-known SFA 
amplitude Q 



Mr, 



dr{^j\r)\H lx {r)\Mr)). (8) 



Here, for times t < r the state of the electron is governed 
by the Hamiltonian Ho , while for t > r its time evolution 
follows the Hamiltonian Hp x . 

The matrix element (|SJ| conveys the following physical 
picture: for times t < r the electron is sufficiently deeply 
bound that to a good approximation its interaction with 
the laser field can be ignored. At time r, it is ionized, and 
the laser intensity is high enough to move the electron so 
rapidly out of the range of the binding potential that now 
the latter can be neglected. 

However, this physical picture is in agreement with the 
formal description only within L gauge. In L gauge, the 
interaction with the laser field is accomplished by the 
scalar potential e$(i) = Hi^t). There is no vector po- 
tential, so that the operator of the velocity is v = p/m. 
Hence, p 2 /(2m) is the operator of the kinetic energy, 
and Ho describes an atom that does not interact with 
the field, even if a field is present. In V gauge, the op- 
erator of the velocity is v = [p — eA(t)]/m where p is 
the operator of the canonical momentum. This is a con- 
served quantity under the dipole approximation, but not 
a physical quantity 1], since p = mv + eA(t), and v is 
a physical quantity while A(t) is not. In consequence, in 
the presence of a laser field, the operator H is not the 
field- free Hamiltonian, and its eigenstate \ipnlto)) does 
incorporate some interaction with the field |8j. Hence, 
in V gauge, the physical picture formulated above is not 
supported by the matrix element JHJl. 

It is instructive to evaluate the matrix element © by 
the method of steepest descent, which is known to work 
very well for sufficiently high intensity. We first recall 
the explicit form of the Volkov wave function 



(r|V$ VJ (*)> 



s -iSp(i) 
(2^ 



%pr 



(x = Y) 



e t(p-eA(t)).r ) ^ = L ) 



(9) 



S p( t ) = ^ f dr(p - eA(r)f 



(10) 



which has the same form in either gauge 9] . 

Via an integration by parts, the transition amplitude 
(jHJ) can be recast in the form 0| 

M p = -i / dr(^)(r)|V(r)|^(r)>, (11) 



which depends on the gauge only via the Volkov wave 
function . Collecting the exponential time dependence 
of the integrand in Eq. 1)11(1 we find that the stationary 
points with respect to r are determined as the solutions 
of the saddle-point equation 



[p - eA(r)] 2 = -2mJ p 



(12) 



The transition amplitude then can be represented as the 
superposition of the contributions of all those solutions 
t s of Eq. I|12|) for which Imi s > 0, with the result 



2ni 



W/ E(i s ) • [p - eA(t s )] 
Only the form factor 



AS p (t s )+I p t s ] 



_ f (p\V(r)\0), (x = V), 



(p-eA(i s )|F(r)|0>, (x = L) 



(13) 
(14) 



depends on the gauge. In V gauge, it is evaluated at 
the momentum p at the detector, which is the same for 
all saddle-point solutions. In L gauge, it is evaluated 
at the instantaneous velocity at the ionization time t s , 
whose component parallel to the laser field according to 
Eq. (|12[) is purely imaginary and can have either sign. 
For a monochromatic linearly polarized laser field, there 
are two solutions t s per cycle of the saddle-point equation 
(|12fl with lmt s > 0, one on either side of the pertinent 
extremum of the vector potential. 

To find out the signs of p — eA(t s ) that correspond to 
the solutions with lmt s > 0, let us consider the vector 
potential A(t) = eAcosujt. We let t = tji + itj, where 
tn and t\ denote the real and the imaginary part of t, 
respectively. The real and imaginary parts of the saddle- 
point equation (|12[> are 



p\\ — eA cos uitn, cosh u)t\ = 0, 



eA sin uit^ sinh u)t\ = \ 2mL 



(15a) 
(15b) 



where pii and p_i_ are the components of p parallel and 
perpendicular to the laser field and the square root may 
have either sign. From Eq. H15al) . the two solutions per 
cycle are such that cos cj<r has the same sign. Then, from 
Eq. I|15b|) and the fact that t\ > for the physical saddle- 
point solutions, we have that sinw% has the opposite 
sign for the two solutions. Hence, the two instantaneous 
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FIG. 1: (color online) SFA electron energy spectrum for emis- 
sion in the (positive) direction of the laser field (four-cycle 
sin 2 -pulse with 4> = 0, lu = 0.056 a.u., Eo = 0.0834 a.u.) in L 
gauge, starting from an initial Is (solid) or 2p (dashed) state. 
The corresponding V gauge result is shown in the inset. 



velocities that enter the L-gauge form factor (|14f> are p — 
eA(t g ) = (±iy / 2mJ p + p^_, pj_). For p^ = 0, they are 
purely imaginary and have opposite sign. This reflects 
the fact that the electric field E(£r) points in opposite 
direction for the two solutions. 

Now, for an even-parity ground state |0), (— a|]/|0) = 
(a|y|0), while for an odd-parity state (— a|y|0) = 
— (a|V|0). Hence, for an odd-parity state, when in 
Eq. IpHfl) the contributions of the two saddle points add in 

V gauge, they substract in L gauge, and vice versa. Con- 
sequently, for an odd-parity initial state, constructive in- 
terference in L gauge implies destructive interference in 

V gauge, and vice versa. In contrast, for an even-parity 
ground state, both gauges predict interference maxima 
and minima at the same positions. As soon as p± ^ 0, 
there is no complete destructive or constructive interfer- 
ence anymore. 

In Figs. ^ 121 an d El we compare the results of the 
SFA in L gauge and V gauge with a numerical solu- 
tion of the TDSE. All calculations have been carried 
out for a 4-cycle linearly polarized laser pulse having 
the intensity 2.4 x 10 14 W/cm 2 (field strength E = 
0.0834 a.u.), and wavelength 800 nm (photon energy u> — 
0.056 a.u.). The electric-field vector is E(t) cos(u>t + <fi)e, 
with the sine-square envelope E(t) = Eq sin 2 j^- for 



< t < T n 



n p T, 



T = 2tt/w, and E{t) = out- 



side this interval. The carrier-envelope phase is (j> = 0. 
Figure n exhibits the results of a numerical computation 
(not using the saddle-point approximation) of the SFA 
amplitude © in L gauge and V gauge, respectively, tak- 
ing fo r \rl>g(t)) the bound state of a zero-range potential 
Hi EH LUl 03 • They illustrate the above statements. In 
other words, in L gauge, everything else being equal, the 
positions of the interference dips for a p ground state 
coincide with those of the interference humps for an s 
ground state. In contrast, in V gauge dips and humps 



FIG. 2: (color online) Same as Fig.0 but computed from the 
numerical solution of the TDSE. 



occur at the same positions regardless of the parity of 
the ground state. Figure |21 presents the corresponding 
TDSE spectrum calculated by methods introduced else- 
where |13J. In order to mimic a short-range potential in 
the TDSE calculations, the Coulomb potential —Z e g/r 
has been cut at r c = 2 a.u. The nuclear charge Z e g was 
adjusted in such a way as to keep the ionization potential 
I p at 0.5 a.u. for both the Is and the 2p state. It has 
been shown recently 14] that the agreement between SFA 
and TDSE low-energy electron spectra improves with de- 
creasing potential range r c . A direct comparison of the 
TDSE and SFA (L gauge) results is presented in Fig. |3] 
The agreement with respect to the energetic positions of 
the various peaks is excellent. Residual discrepancies are 
observed in the shape of the spectrum for low energies, 
especially for the p ground state, and are likely due to 
the different large- distance behavior of the wave func- 
tions (zero range for the SFA vs. cut Coulomb for the 
TDSE). 

The exact solution for ionization of negatively charged 
ions that is available in the context of effective-range 
theory exhibits the interference dips in complementary 
positions for s and p ground states [l^, in agreement 
with the L-gauge SFA and the exact TDSE solution. 
The L-gauge SFA also appears to be supported by the 
experimental data: the above-threshold-detachment en- 
ergy spectrum for the negative F~ ion |lfi| |. which has 
a p ground state, displays a pronounced change of its 
slope at the energy where the L-gauge SFA predicts an 
interference dip [ToL ITT| . 

For elliptical polarization, for ellipticities higher than a 
certain critical value the saddle-point equation (|12fl only 
has one solution per cycle rather than two, so that the 
interference ceases to exist [T^ |. This is so, in particular, 
for circular polarization. Recently, the latter case was 
considered in detail |l8j]. Even in the absence of inter- 
ference, the form factor (|14f> is still different in L gauge 
and in V gauge. For an s ground state |0), the form 
factor {p| V|0) has a maximum for p = and decreases 
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FIG. 3: (color online) Direct comparison of TDSE and SFA 
(L gauge) photoelectron spectra. The initial states are (a) Is 
and (b) 2p, the other parameters as in Fig.Q For convenience, 
the TDSE spectra are rescaled (but unshifted in energy). 



with increasing |p|, while for a p state, it has a zero at 
p = and extrema away from p = 0. In Ref. [l8| . 
for ionization of F~ by a circularly polarized laser field, 
the energy spectrum was calculated in either gauge. The 
V-gauge spectrum peaks at a higher energy than the L- 
gauge spectrum, which conforms with the considerations 
given above. Moreover, Wigner's threshold law is only 
reproduced in L gauge 

Before concluding, we recall that in a numerical solu- 
tion of the TDSE the choice of gauge is "merely" a ques- 
tion of convenience. Generally, convergence is faster in 
V gauge where fewer angular momenta contribute, much 



faster indeed for high intensity and low frequency p^j . In 
contrast, in approximations such as the SFA, the choice 
of gauge is a contributing factor for the quality of the 
approximation. In fact, making formally the same ap- 
proximation in two gauges may correspond to different 
approximations physically. A general argument in favor 
of the Lgauge for use in the SFA has been put forward 
in Ref. |6j : the L-gauge interaction Hamiltonian @ puts 
the emphasis on large distances from the atom, where 
the Volkov wave function is a good approximation to the 
final state. In a similar vein, we add that it appears to 
make more sense to evaluate the form factor (|14|) at the 
instantaneous velocity at the time of ionization (as in L 
gauge) rather than at the drift velocity (as in V gauge), 
which for low frequencies the electron does not assume 
before it is far away from the ion. 

On the basis of a comparison with the solution of 
the time-dependent Schrodinger equation, we conclude 
that the strong-field approximation applied to above- 
threshold detachment of negative ions affords a better de- 
scription in length gauge than in velocity gauge. In view 
of the fundamental significance of the SFA for strong-field 
physics, it is of great importance to find out which gauge 
is better suited for above-threshold ionization of atoms 
and molecules as well as nonsequential double ionization. 
In all of these cases, the two gauges are known to yield 
different answers as well. 

We enjoyed discussions with M. Yu. Ivanov, H. 
G. Muller, and H. R. Reiss. This work was sup- 
ported by VolkswagenStiftung, Deutsche Forschungsge- 
meinschaft, and the NSERC Special Opportunity Pro- 
gram of Canada. 
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